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Abstract 

The one and two-particle form factors of the energy operator in the two-dimensional 
Ising model in a magnetic field at T = T c are exactly computed within the form 
factor bootstrap approach. Together with the matrix elements of the magnetisation 
operator already computed in ref. 0, they are used to write down the large distance 
expansion for the correlators of the two relevant fields of the model. 



The last years have seen important progresses in the non-perturbative study of two- 
dimensional quantum field theories and related statistical mechanical models. If confor- 
mal symmetry provided us with an exact description of critical points and universality 
classes |2|], the study of off-critical models turned out to be better approached within 
the framework of relativistic scattering theory. In fact, if the off-critical model under con- 
sideration is integrable (i.e. admits infinite conservation laws), it can be usually solved 
exploiting very general bootstrap techniques |3], f|, ||, |6| . This circumstance appears to be 
particularly important in light of the fact that a large number of physically interesting 
two-dimensional systems can actually be described in terms of integrable models. A re- 
markable example is provided by the scaling limit of the two-dimensional Ising model in 
a magnetic field at T = T c (IMMF in the sequel) 0] . It can be formally described by the 
action 

A = Ac ft + hjd 2 x a(x) , (1) 

where Ac ft denotes the action of the conformal minimal model Ai^^ and a(x) the mag- 
netisation operator of scaling dimension 2A°" = 1/8. The coupling constant h (magnetic 
field) has physical dimension h ~ m 15 / 8 , m being a mass scale. Apart form the mag- 
netisation operator, the only other relevant scaling field in the Ising model is the energy 
density e(x) with scaling dimension 2A e = 1. 

Zamolodchikov showed that the theory described by the action ([!]) possesses an infinite 
number of integrals of motion which can be used in order to determine the exact particle 
spectrum and S'-matrix of the theory ||]. He found that the spectrum consists of eight 
massive particles A a (a = 1, 2, . . . , 8) whose masses stay in the following ratios with the 
mass mi of the lightest particle 
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The interaction between these particles is described by a factorised, reflectionless S'-matrix 
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characterised by the two-particle amplitudes 

tanh | (9 + ina) 1 . , 

(3) 



s ab (9)= n 
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tanh | [6 — irca) 



The set of numbers A a b and the multiplicity factors \i a are given in Table 1. We use the 
standard rapidity parameterisation of the on-shell momenta p% = (m a cosh 9 a , m a sinh 6 a ), 
so that 6 = 6 a - 9 b in (|). 

In ref. the knowledge of the S-matrix (|]) was exploited in order to approach the 
computation of the correlation functions of the model (ED) within the form factor boot- 
strap method. The basic idea of this approach is to express the (euclidean) correlation 
functions (e.g. the two-point ones) as a spectral sum over a complete set of intermediate 
multiparticle states 



G^ 2 {x) = ($i(z)$ 2 (0)) 
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x (0|* 1 (0)|A ll (fl 1 )...^(e»)><An(^---^.(fi»)|$ 2 (0)|0> , 
and to exploit the fact that the form factors (FF) 

Ft.. an {0i, ■ • • A) = (0\*(0)\AaM ■ ■ ■ AaM) ■ (5) 

are exactly computable in integrable models once the S"-matrix is known || ^ . 

The spectral series @ is manifestly a large distance expansion. Nevertheless, it has 
been observed in several models [9-13] that it is characterised by a fast rate of convergence 



also at intermediate and short distance scales (see in particular ref. [ID] for a theoretical 
justification of this property), so that a truncation of the series including few lowest states 
turns out to be sufficient for most practical purposes. In ref. the one and two-particle 
FF of the magnetisation operator in the IMMF were exactly computed and shown to be 
sufficient to reproduce with remarkable accuracy the numerical data for G aa (x) avalaible 



from Monte Carlo simulations |TJJ. It is the purpose of this letter to carry out a similar 



program for the FF of the energy operator e(x) and to write down the large distance 
expansion for the correlators G e£ (x) and G ue (x). 

Let's briefly recall the basic strategy for the computation of FF in the IMMF referring 
the reader to ref. for details. Form factors can be generally computed in an integrable 
model using a recursive procedure based on a set of residue equations relating matrix 
elements with different particle content. For instance, if the scattering amplitude S a b(9) 
has a simple pole with positive residue (r^ fc ) 2 at 6 = iu c ab corresponding to the particle 
A c with mass ml = m 2 a + m 2 + 2m a m b cos u° ab appearing as a bound state in the direct ab 
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channel, then we can write 

F*{9 * iul h ) * 1 ^f-F? , (6) 

as well as similar relations among higher matrix elements containing spectator particles. 
Other recursive equations are associated to the higher poles in the S-matrix and to 
the "kinematical" poles in the matrix elements. Together with the equations ruling the 
monodromy properties of FF, these residue equations provide a system of linear relations 
whose general solution amounts to a complete classification of the operator content of the 
theory |15] . 

In theories with diagonal scattering, the basic information about the structure of FF is 
already encoded in the two-particle matrix elements. In the IMMF they can be generally 
parameterised as |7j [] 

QUO) 



FM = ^fj^F2 in {9) , (7) 



where 
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The terms F^ in (8) and D a b(8) in eq. (0) take into account the monodromy properties 
and the singularity structure of the matrix elements, respectively, and are both uniquely 
determined by the knowledge of the 5-matrix. The whole information about the operator 
is then contained in the polynomial 

Nab 

Qa6(0) = £4cosh fc 0, (13) 
k=0 

whose (operator dependent) coefficients c k ab are the only remaining unknowns. 
1 In this letter we consider only scalar operators. 
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The identification of specific operators out of the general solution of the recursive 
equations is a nontrivial task, especially in theories lacking any internal symmetry, as is 
the case for the IMMF. A progress with respect to this problem was made in ref. 0] where 
it was shown that in a unitary theory the FF of an operator satisfy the asymptotic 
bound 

lim (9i, . . . A) < const. e A *^ , (14) 

\0i\^oo 

2A* being the scaling dimension. The consequences of this result are easily illustrated 
using the parameterisation (|7]). Consider the simplest two-particle FF in the IMMF, 
F*(6>). Since G a {6) ~ exp(|0|/2) as \6\ — > oo, the constraint fll4]) implies that, for any 
relevant scalar field f(x) (A v < 1), the total degree Nu of the polynomial Qfi(0) must 
be less than 2. Since it can be checked that no solution of the residue equations exists if 
N n = 0, one concludes 

Qii(0) = c l n coshe + c° n , 4^0. (15) 

It turns out that, once the FF Ff^iO) corresponding to a specific relevant operator <p(x) 
has been fixed assigning the two coefficients c\ ± and c° 1( all the FF of <p(x) can be uniquely 
determined using the residue equations. This amounts to say that the solutions of the FF 
bootstrap for the relevant scalar operators of the IMMF form a two-dimensional linear 
space, which is what expected from the fact that such operators can only correspond 
to linear combinations of a(x) and e(x). Solutions corresponding to operators not just 
differing for an inessential normalisation constant can be labelled by the ratio 

z^ C \. (16) 
c n 

It is of particular physical interest to determine the values of this ratio which select the two 
scaling fields a{x) and e{x). The problem for a{x) was solved in ref. [|7| using the relation 
between the field perturbing the conformal point and the trace of the energy-momentum 
tensor 

Q(x) = 27ih(2-2A a )a(x) , (17) 

and exploiting the constraints imposed on the matrix elements of Q(x) by energy-momentum 
conservation. On the other hand, no similar method can be used for e(x) and an alterna- 
tive way must be found in order to characterise the FF of this operator. 

The fields cr(x) and e(x) are uniquely identified by the short distance behaviour of 
their correlation functions predicted by the conformal operator product expansion |3] 

x\ -> 

^0 (18) 
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Since in the FF approach the determination of the exact ultraviolet behaviour of the 
correlators requires in principle the resummation of the spectral series (^), it seems quite 
difficult to make a direct use of eqs. (|18D for the determination of z a and z £ . Nevertheless, 
eqs. fll8|) suggest that, if a property characterising the FF of the scaling fields a(x) and 
e(x) exists, it is probably related to the high energy asymptotics of the matrix elements. 
Moreover, if such property should be able to select the scaling fields among their linear 
combinations, it must be non-linear in the operator. Interestingly enough, a property 
with these features is known in the FF literature. It was noticed in the context of the 
sine-Gordon model that the FF of the operators exponential of the elementary field satisfy 
the cluster property || |TB|, |T7| 



Ft..a^ ... an • • • , Ok+a, 9 k+x , . . . A) = K..*> (*i> • • • > **)*2L „.«. Vh+i, ■ ■ ■ A) 

(19) 
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where 



Ft..aS^ -,e n ) = ^Fl^, ...,e n ). (20) 

It is known that massive deformations of minimal models can be obtained from sine- 
Gordon through a suitable restriction of the Hilbert space and that some exponential 
operators are mapped into the scaling fields of the restricted models. In refs. [16|, [T3], 
it was shown for some specific cases that the factorisation property flO^D survives the 
reduction procedure and is satisfied by the FF of the scaling fields in the reduced models. 
Here we simply assume that the cluster property flTP| ) characterises the FF of the operators 
o~(x) and e(x) in the IMMF and provide what we think is strong evidence that this is 
indeed the case. 

Before proceeding further, notice that the asymptotic relation ( |T9"D involves the vev 
($). In the general case, the computation of this quantity is a nontrivial problem (the 
thermodynamic Bethe ansatz only provides the vev of the field which perturbs the confor- 
mal point [[BJ). It is then remarkable that, if the FF of $(x) satisfy the cluster property, 
the vev can be obtained, for instance, as 

Going back to the determination of z a and z £ , consider the two-particle FF Fi 2 (8). Follow- 
ing the same arguments used above for Ff^O) one concludes that Qi 2 (0) is a polynomial 
of degree 2 in cosh^. Since the amplitudes Sn(9) and Si 2 (9) have common poles corre- 
sponding to the particles Ax, A 2 and A 3 (see Table 1), eq. @ provides the linear system 

l-Rese^F^O) = ^rRes e=iu c 2 F? 2 (e) , c = 1, 2, 3 . (22) 

1 11 1 12 
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Once an overall normalisation of the operator <f(x) has been fixed, these equations 
uniquely determine the coefficients c\ 2l c\ 2 and c\ 2 in terms of z v . Finally, in order 
to search for solutions satisfying the cluster property, we use eq. ([H]) and require 

Fl F2 (23) 



lim^F^) Mmg^F^e) ' 

with Fi and F 2 also determined in function of z v using eq. (||). 

There exist only two values of the parameter z v satisfying the last equation. One of 
them exactly coincides with the value of z a which had been determined in ref. [0] without 
any reference to the cluster property 

, CT = ^i±^ = 4.869840.. . (24) 
2mf 

As we said above, once this initial condition has been fixed, all the matrix elements of 
a[x) can in principle be determined using the residue equations only, without further use 
of the cluster property. All the one-particle and several two-particle FF of a(x) are given 
in tables 2 and 4, respectively 0. It must be stressed that all the two-particle FF of a(x) 
computed in this way automatically satisfy the cluster property which then should be 
regarded as characteristic of the whole solution selected by the initial condition ( f2"4|) 0. 
The same pattern is observed for the solution arising from the other value of z v satisfying 
eq. (|23|), which we identify with z £ f] 

z £ = 1.255585.. . (25) 

The one and two-particle FF corresponding to this initial condition are contained in tables 
3 and 5; they were used in ref. |2(J in order to compute the corrections the energy spectrum 
of the theory (H) undergoes under a small thermal perturbation induced by the energy 
operator e(x). The remarkable agreement observed in ref. []2D| between the theoretical 
predictions and the data coming from a numerical diagonalisation of the Hamiltonian 
strongly supports the conclusion that the cluster property correctly selects the matrix 
elements of both relevant scaling operators in the IMMF. 

The results contained in tables 2-5 can be used in the spectral representation (Q) in 
order to write down the large distance expansion of the correlators G aa (x), G ae (x) and 

2 Of course, these results coincide with those of ref. 0. 

3 In other words, the ratio in eq. (^l|) is the same for any a and b. Moreover, once the standard CFT 
normalisation of the operator is adopted, the value obtained for (a) exactly coincides with the result 
provided by the thermodynamic Bethe ansatz fl9|| . 

4 The analytic expression for z E we dispose at the moment looks complicated and uninspiring and we 
prefer to quote only its numerical value. 
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G ee (x). The leading infrared contributions are simply 

G^ 2 (x) = ($ 1 )($ 2 ) + - £ F^F**K (m a \x\) + O (e^M) (26) 



7r o=l 



where Kq(x) is the modified Bessel function. In ref. 0, the fast convergence of the FF 
series for G aa {x) was tested against the Monte Carlo data available for that correlator; 
as far as we know, no similar data exist for G ae (x) and G ££ (x) but it seems natural to 
expect a similar convergence pattern. An integral check is provided by the following sum 
rule for the scaling dimension pi] 

A* = --1-r / d 2 x(Q(x)m)c ■ (27) 



4tt($) 

The results obtained for A a and A £ using in the last formula the spectral representa- 
tion of the correlators G aa (x) and G ae (x) are contained in tables 6 and 7, respectively 
(A* an denotes the contribution coming from the intermediate state containing the par- 
ticles A ai . . . A an ). We recall that the exact results are A°" = 0.0625 and A e = 0.5. The 
reason for the slower convergence of the series for A e appears quite clear. Indeed, accord- 
ing to (|T8D, G aE {x) is much more singular than G aa (x) as x — > 0. As a consequence, the 
integral in eq. ( p7|) receives a larger contribution from short distances for $ = e than for 
$ = a. On the other hand more and more terms in the spectral series are nedeed in order 
to approximate precisely the correlators at small x. 

In conclusion, it is clear that it would be higly desirable to reach a satisfactory physical 
understanding of one of the basic ingredients we used in this letter, namely the cluster 
property (|19|). It is very tempting to argue that this property (or better, a suitable 
generalisation applying also to theories with internal symmetries) characterises the matrix 
elements of the scaling fields in two-dimensional quantum field theories. The Smirnov's 
observation that, if the FF of an operator factorise asymptotically as in (|T5D, then it 
is particularly simple to show that the correlator G$$(x) behaves as a power law at short 



distances [16], seems to go in this direction. We hope that the results of this letter will 



stimulate further investigations on this point. 
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Table Caption 



Table 1 . Two-particle scattering amplitudes of the IMMF. Each factor (7)^ stays for 
tanh \ {p + ^^0) / tanh \ {0 — m M . The indeces i placed above the functions 
(7) correspond to the particles Ai appearing as bound states in the ab channel. 

Table 2 . One-particle Form Factors of the operator a(x). The results are given in units 
of m}^ 8 and refer to the normalisation of the operator in which (a) = m{^ 8 . 

Table 3 . One-particle Form Factors of the operator e(x). The results are given in units 
of m\ and refer to the normalisation of the operator in which (e) = m\. 

Table 4 . Coefficients of the polynomials Q a ab (6). The results are given in units of m}^ 8 
and refer to the normalisation of the operator in which (a) = m^ 8 . 

Table 5 . Coefficients of the polynomials Q £ a b(@)- The results are given in units of m\ 
and refer to the normalisation of the operator in which (e) = m\. 

Table 6 . The first eight contributions to the sum rule for A°\ 

Table 7 . The first eight contributions to the sum rule for A e . 
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c 15 — 


-260.7643072 


r 3 - 
c 15 — 


-4719.877128 


r 2 - 

c 15 — 


-15172.07643 


c 15 — 


-17428.22924 


c 15 — 


-6716.787925 


r 4 - 
c 23 ~~ 


-92.73452350 


r 3 - 
c 23 ~~ 


-1846.579035 


r 2 - 
c 23 ~~ 


-6618.297073 


c 23 ~~ 
c 23 ~~ 


-8436.850082 
-3579.556465 


c 33 ~~ 


-1197.056497 


r 4 - 
c 33 — 


-30166.99117 


r 3 - 
c 33 — 


-150512.4122 


r 2 - 
c 33 — 


-301093.9432 


c 33 ~~ 


-267341.1276 


c 33 — 


-87821.70785 



Continued 
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r 6 — 
c 25 ~~ 


1425.995027 


c 25 ~~ 


44219.03877 


r 4 - 
c 25 ~~ 


286184.1535 


r 3 - 

c 25 ~~ 


788413.2178 


r 2 - 

c 25 ~~ 


1078996.488 


c 25 ~~ 


725356.4417 


c 25 — 


191383.5734 


r 5 - 
c 17 — 


190.8548023 


r 4 - 
c 17 — 


4633.706068 


r 3 - 
c 17 ~~ 


21406.72691 


c 17 — 


39514.82959 


c 17 = 


32456.91939 


c 17 ~~ 


9906.265607 




-7249.785565 




-276406.7236 


r 5 - 

c 44 ~~ 


-2299573.212 


r 4 - 

c 44 _ 


-849276.3526 


c 44 _ 


-16615618.39 


r 2 - 

c 44 _ 


-17950817.11 


C44 = 


-10139089.36 


c 44 — 


-2341590.241 
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-70.00917205 


c°n = 


-87.90247670 


r 2 - 

c 12 ~~ 


-466.3008246 


C 12 = 


-1307.331521 


c 12 — 


-853.2803886 


r 3 - 
c 13 ~~ 


-43021.45153 


r 2 - 

c 13 ~~ 


-182413.2733 


c 13 = 


-241929.7678 


c 13 ~~ 


-102574.1349 


r 3 - 

c 22 — 


-2193.896354 


r 2 - 

c 22 — 


-10870.05277 


C 22 = 


-16161.44508 


c 22 — 


-7510.235388 




= 2074.636471 


r 2 - 


= 9881.413381 




= 14357.04570 




= 6568.762583 


r 4 - 
c 15 — 


-30333.56619 


r 3 - 
c 15 — 


-198757.2340 


r 2 - 

c 15 — 


-447504.5720 


c 15 — 


-422808.9295 


c 15 — 


-143743.2050 


r 4 - 
c 23 ~~ 


-11971.94909 


r 3 - 
c 23 ~~ 


-81253.72269 


r 2 - 
c 23 ~~ 


-186593.8661 


c 23 ~~ 
c 23 ~~ 


-178494.3378 
-61194.62416 


c 33 ~~ 


-195385.7662 


r 4 - 
c 33 — 


-1743171.802 


r 3 - 
c 33 — 


-5603957.324 


r 2 - 
c 33 — 


-8422606.859 


c 33 ~~ 


-6035102.896 


c 33 — 


-1668721.004 
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r 6 - 

c 25 ~~ 


289831.4882 


r 5 - 

c 25 ~~ 


3275586.983 


r 4 - 
c 25 ~~ 


13872077.63 


r 3 - 

c 25 ~~ 


29236961.96 


r 2 - 

c 25 ~~ 


32979257.31 


c 25 ~~ 


19100224.04 


c 25 ~~ 


4471623.121 


r 5 - 
c 17 — 


30394.23374 


c 17 — 


274294.8033 


r 3 - 
c 17 ~~ 


897781.3229 


c 17 — 


1.375919456 


c 17 = 


1.004969466 


c 17 ~~ 


282938.1974 




-1830120.693 




-25699492.93 


r 5 - 

c 44 ~~ 


-138411873.8 


r 4 - 

c 44 _ 


-384776478.8 


c 44 _ 


-608371427.1 


r 2 - 

c 44 _ 


-553818699.0 


C44 = 


-270964337.7 


c 44 — 


-55283137.91 
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A? 


0.0507107 


AS 


0.0054088 


A CT 


0.0010868 


An 


0.0025274 


K 


0.0004351 


A? 2 


0.0010446 


A^ 


0.0000514 


^13 


0.0002283 


Ao- 
partial 


0.0614934 



Table 6 



Af 


0.2932796 


A! 


0.0546562 


A e 

^3 


0.0138858 




0.0425125 


A £ 

4 


0.0069134 


Af 2 


0.0245129 


Af 


0.0010340 


A £ 

^13 


0.0065067 


A e 

partial 


0.4433015 
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